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ABSTRACT 


Interpolating  elastica  are  the  extremals  for  the  functional  j c (s)ds,  which 

0 

is  the  integral  of  the  square  of  the  curvature  with  respect  to  arc  length,  in  the  feunily 


of  plane  curves  that  interpolate  at  (not  prescribed)  arc  lengths  ®o  ^ ^1 


< s 


a^rescribed  configuration  of  points  P„»P. 


If  at  one  or  both  terminals  the 


indefinite  length  with  sleeve  supports  anchored  at  P_,P, ,...p  , 

0 1 n 


■0'^l'’“'^n' 

slope  is  prescribed,  the  extremal  is  said  to  be  angle-constrained,  otherwise  free.  The 
curvature  functional  represents  the  elastic  strain  energy  of  a thin  elastic  beam  of 

which  allow  the  beam 

to  slide  through  without  friction  and  to  rotate  freely  (except  at  the  end  supports  if 
angle-constrained) . The  interpolating  elastica  are  also  known  as  nonlinear  spline 
curves.  It  is  known  that  the  infimum  of  the  strain  energy  is  0 in  all  cases,  hence 
cannot  be  attained  if  the  points  Pr,»P, »-..,P  do  not  lie  on  a ray.  On  the  other  hand, 
interpolating  elastica  are  known  to  exist  for  a variety  of  configurations,  and  this 
report  investigates  whether  these  extremals  make  the  strain  energy  a local  minimum  or 
not  (i.e.  vihether  they  are  "stcdile"  or  "unstcible") . Several  general  stcd^ility  criteria 
are  established  and  they  are  used  to  decide  the  stability  of  some  specific  elastica. 
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SIGNIFICANCE  AND  EXPLANATION 


It  is  an  old  technique  of  draftsmen  to  use  a mechanical  spline  to 
pass  a smooth  curve  through  a prescribed  set  of  points  in  a plane.  Curves 
which  are  obtained  in  this  way  (interpolating  elastica,  also  called  non- 
linear spline  curves)  may  be  considered  as  the  equilibrium  positions  of 
thin  elastic  beams  which  are  constrained  to  pass  through  short,  friction- 
less, freely  rotating  sleeve  supports,  anchored  at  the  interpolation 
points.  The  strain  energy  of  such  a beam  is  given  by  the  integral  of 
the  square  of  the  curvature  with  respect  to  arc  length,  and  equilibrium 
requires  that  the  position  be  such  that  the  energy  be  minimal  for  the 
given  interpolation  conditions.  However,  a global  minimum  cannot  be 
attained  (except  in  the  trivial  case  of  the  unbent  beeim)  since  the  energy 
can  be  made  arbitrarily  small  by  using  sufficiently  large  loops  between 
the  supports.  Instead  one  looks  for  local  minima  which  guarantee  sta- 
bility against  small  perturbations.  In  this  report  some  general  stability 
criteria  are  established  and  some  specific  interpolating  elastica  are 
investigated  for  stability.  Except  for  a few  previous  isolated  observa- 
tions these  seem  to  be  the  first  proven  results  on  the  stability  of  inter- 
polating elastica. 


\ 


The  responsibility  for  the  wording  and  views  expressed  in  this 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


STAIUI>ITY  OF  INTKHPOLATINr,  KI^STICA 


Mi  char*  1 Oolomi^) 


1 . 2f*l-  ion 

Klastica  aro  tho  piano  curves  with  "normal  roprosontation"  s 0(.s)  (s  denoLos 
arc;  Ior»«jth  and  0(s)  the  am/Io  of  i nc*  1 i nat ic^n  at  s)  which  are  solution.s  of  fho 
di  f ferontiai  oc^uation 

(l.l)  ~<)'^(S)  = A(sin(0  - Oj)  - u| 

where  A,0^  and  u arc  real  constants  (see,  e.cj.  H,  Article  26^1).  (1.1)  is  the 

Fuler  e<piation  lor  the*  variaticjnaJ  problem 

s s s 

(K2)  'S  J O'  “ 0,  J cosO  1=  b,  / sinO  = d 

0 0 0 

wliere  s,b,<i  are  prcjscribed  (sec.*  above  referemee  or  (2,  Proj^.  1.2)).  'n»e  intc^cjral 

s 2 

/ O'  ' represents  (with  the  f>rc}per  choice  of  units)  the  strain  onerqy  of  a thin 

0 

elastic  bc*am  of  uniform  cre^ss  section  of  lenqth  s,  and  the  side  conditions  in  (1.2) 
spc:ci  fy  the  relative;  f>osition  of  the*  ends  of  the  bent  beam. 

'file  elastica  described  by  (1.1),  when  considered  for  all  values  of  s,  have 
infinitely  many  inflection  (xjints,  0*  (s)  = 0 when  sin(n(s)  - 0^)  = n,  and  are 

therefore  called  inflectional  elastica  (see  (1,  loc  cit.l).  Below  we  will  consider 
only  elastica  for  which  a » 0;  geometrically  sf^ealcinq,  these  are  curves  for  which 
the  variation  of  0 betwc.*on  consecutive  inflection  jxiints  is  ti.  We  refer  to  them 
as  simp  If.*  elastica.  All  the  simjile  elastica  are  obtained  from  a particular  one  by 
similarity  transformations. 

Tlic  inter£>olatinq  elastica  (so  named  by  M.  A.  Malcolm  in  (11)  cernsist  of  finitely 

nwiny  subarcs  of  the  simple  elastica,  fitted  toejether  so  that  a smooth  curve  with 

c.ontlnuous  curvature  results  which  has  jump  discontinuities  of  the  derivative  of 

the  curvature  only  at  the  "knots"  Such  an  interpolatinq  elastica  H 

1 n-l 
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with  normal  representation  9 is  the  solution  of  the  variational  problem 


I 

1 


t 


; 

i 


I 


(1.3) 


S / 0'^  =0,  / cos6  = b.,  / sinB  = d (i  = 

' 1 1 


s . 

1 


®i-l  ®i-l 


whore  the  are  prescribed  (they  are  the  coordinates  of  the  vector  p^  - ^), 

but  the  arc  lengths  ^ . . . < s^  of  the  terminals  Pq<P^  of  the  knots 


. ,p  are  varied  (see  (2,  loc.  cit.]  or  (3,  Sec.  21).  If  the  ends  P-.P 


1 ‘^n-l 


0'*^1 


are  "free"  then  the  natural  boundary  conditions 

(1. 3a)  0'  (0)  = 0,  0'  (s  ) = 0 

n 

are  appended  to  (1.3).  Frequently  we  shall  be  concerned  with  "anqle-constrained" 

interixjlating  elastica;  in  this  case  we  are  given 

(1. 3b)  0(0)  = a,  0(s  > = B . 

n 

The  solutions  of  (1.3),  (1.3a)  represent  possible  equilibrium  positions  (stable  or 
not)  of  a thin  elastic  beam  of  indefinite  length  which  is  constrained  to  pass  through 
frictionless  freely  rotating  small  sleeves  anchored  at  the  positions  Pg'^i' ' ‘ ' 'Pn' 

If  the  sleeves  at  the  terminals  Pg'P„  ^te  pinned  then  (1.3b)  replaces  (1.3a).  The 

interfralating  elastica  are  a reasonable  mathematical  model  for  the  mechanical  spline 
used  by  draftsmen  to  pass  a smooth  curve  through  the  given  points  p^.p^^, . . . ,p^. 

They  are  also  called  nonlinear  (interpolating)  splines  (see,  e.g.,  (4))  and  were 

referred  to  as  extremal  interpolants  for  the  configuration  (p_,p ,P  } in  [21. 

0 1 n 

We  still  will  refer  to  them  by  this  name  in  the  sequel. 

The  solutions  of  (1  3)  are  definitely  not  absolute  minima,  except  in  the  trivial 

case  where  P„,P,,...,P  lie  (in  this  order)  along  a ray  (and  moreover,  a » B = 0 
0 1 n 

in  case  of  end  conditions  (1.3b)).  This  was  first  pointed  out  by  the  authors  of  (5). 
s 

n 2 

/ 0'  can  be  made  arbitrarily  small  by  using  large  interpolating  circular  loops. 

0 
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The  sulutioiia  are  often  toforreti  to  as  local  minima,  altliouqh  no  proofs  are  qiven 
that  they  are  indeed  extrema  of  this  kind.  Only  in  (2,  Theorem  6.1)  was  it  proved 
that  tiie  nontrivial  simple  elastica  i nterpolatinq  2 points  are  nonstable,  i.e.,  they 

®1  2 

do  not  reini'sent  local  minima  of  / O'  . It  is  the  objective  of  this  paper  to 

establish,  for  several  known  extremal  interpolants , whether  they  are  local  minima 
or  not  (stable  or  unstable). 

The  fact  that  the  extremal  intcrpolants  do  not  represent  minima  nor,  in  qeneral, 
local  minima  of  the  functional  / 0*^  is,  probably,  the  major  reason  for  the  lack 
of  qen€*ral  existence  results  and  of  qood  computational  procedures.  (For  an  existence 
proof  limited  to  lenqth-restricted  extremals,  see  [6).  In  [7,  Theorem  3)  it  is 
proved  that  if  there  is  a length- restricted  extremal  of  ’’unstable  length",  there  is 
also  an  interpolating  local  minimum,  but  no  nontrivial  length-restricted  extremal  of 
unstable  length  is  presented.  Existence  of  length-prescribed  extremals  and  of 
iinre.^itrictcc)  extremal  intcrpolants  close  to  a ray  intorpolant  is  proved  in  [2,  Appendix 
and  Theorem  7.4),  where  also  many  examples  of  specific  interpolants  are  given,  which 
were  not  known  before.  For  a survey  of  old  and  new  computational  procedures,  see  (3).) 
In  the  discussion  of  stability  (that  is,  whether  the  extremals  are  local  minima  or  not) 
we  naturally  restrict  ourselves  to  cases  wherj  existence  of  extremal  interpolants  has 
been  proved  or  is  postulated. 

In  Section  2 the  variational  equations  for  interixalating  splines  in  normal 
representation  are  derived,  without  recourse  to  Lagrange  multiplier  theory,  and  as 
a preparation  for  the  computation  of  the  second  variation.  In  Section  3 the  second 
variation  is  used  for  stability  criteria  (Jacobi*s  condition):  an  explicitly  given 
quadratic  functional  must  be  positive-definite,  or  equivalently,  a nonconventional 
linear  second-order  boundary  value  problem  must  have  only  positive  eigenvalues.  In 
Section  4 it  is  proved  that  interpolating  splines  close  (in  a precise  sense)  to 
stable  ones  are  stable  and  those  close  to  strongly  unstable  ones  are  unstable.  This 
result  is  then  used  to  prove  that  splines  that  interpolate  configurations  close  to 


a ray  configuration  (whose  existence  was  proved  in  (21)  are  stable  (even  in  the  case 
of  free  terminals) . This  is  probably  the  first  general  existence  proof  for  locally 
minimizing  interpolants  which  are  not  length-restricted.  In  Section  h it  is  proved 
that  the  extremal  2-point  interpolant  consisting  of  n ^ 1 complete  loops  of  the 
simple  elastica  is  unstable  even  if  angle-constrained  (in  [27!  the  instability  was 
proved  for  the  free  elastica) . If  the  angle-constrained  2-point  interpolants  is  a 
proper  subarc  of  one  loop  of  the  simple  elastica  (hence  has  no  inflection  point)  then 
it  is  stable,  and  any  angle-constrained  2-point  interpolant  that  contains  one  complete 
loop  of  the  simple  elastica  is  unstable.  The  proof  for  these  last  results  is  contained 
in  Section  6;  it  is  built  mainly  on  the  discovery  of  the  eigenfunction  iielonging  to 
the  eigenvalue  0 for  the  second  variational  equation  that  goes  with  the  one-loop 
angle-constrained  simple  elastica.  By  an  extension  of  this  method  it  is  proved  in 
Section  7 that  if  an  angle-constrained  interpolant  contains  an  interior  inflection 
fjoint  then  it  is  stable  if  it  contains  neither  the  left  nor  the  right  "stability 
focus".  These  are  points  on  the  simple  elastica  which  are  situated  symmetrically  with 
respect  to  the  inflection  point,  not  far  from  the  neighboring  inflection  points.  If 
the  angle-constrained  2-{joint  interpolant  with  one  inflection  point  contains  both 
stability  foci  it  is  unstable.  The  general  result  on  the  stability  of  such  2-point 
interpolants  is  stated  with  the  use  of  what  we  call  "conjugate  points".  If  p is 
a point  on  a simple  elastica  arc  containing  one  inflection  point  there  is  a conjugate 
point  defined  by  a transcendental  equation,  and  it  is  also  given  a geometric 

interpretation  (p  and  p^^  are  on  opposite  sides  of  the  inflection  point;  if  p is 
a stability  focus  then  p^  is  the  other  stability  focus) . The  angle-constrained 
elastica  is  stable  if  and  only  if  it  contains  no  pair  of  conjugate  points.  If  the 
2-point  extremal  interpolant  is  free  at  one  end  and  angle-constrained  at  the  other 
end,  then  it  is  stable  if  and  only  if  it  contains  no  stability  focus.  Section  8 
contains  the  most  important  stability  results.  It  is  first  proved  that  a necessary 
condition  for  the  steibility  of  extremal  N-point  interpolants  is  that  each  arc  between 
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consecutive  nodes  be  "proper",  i.e.  internal  arcs  do  not  contain  a pair  of  conjugate 
points,  and  the  terminal  arcs  do  not  contain  a stability  focus.  Then  a computable 
"stability  function"  of  (N-2)  real  varieibles  is  defined  for  the  extremal  N-point 
interpolant  under  investigation,  which  has  a critical  value  at  the  point  that 
corresfxsnds  to  the  extremal.  It  is  proved  that  the  extremal  is  stable  if  and  only 
if  the  critical  value  is  a local  minimum.  These  results  are  applied  to  decide  the 
stability  of  some  3-point  and  4-point  extremal  interpolants . In  this  connection  it 
is  also  shown  that  the  often  repeated  claim  (first  appearing  in  15))  that  a certain 
4-point  configuration  has  no  interpolating  elastica  is  false.  In  the  last  section 
we  show  that  the  closed  extremals  which  interpolate  the  verticc"  of  a regular  n-gon 
(n  / 3)  (their  existence  is  proved  in  [2,  Sec.  81)  are  stable. 


Let  s 


er-Lagrange  conditions  for  the  interpolating  spline  in  normal  representation. 

'<•5(0),  O^s^s  be  the  normal  representation  of  an  admissable  inter- 

polant  C for  the  configuration  {p„/P, }•  Here  s denotes  the  arc  length 

0 1 n 

along  the  curve  and  0(s)  the  angle  that  C malces  at  arc  length  s with  a 

reference  line.  The  interpolation  conditions  are 


/ ^ cos8(s)ds  = b^,  / sine(s)ds  = d^,  i = l,...,n 


where  b.,d.  are  given  numbers,  and  the  nodes  0 = s_  < s,  < . . . < s = s are  the 
11  0 1 n 

arc  lengths  at  which  C passes  through  the  interpolation  points  Pg'^l' ' ‘ ' '^^n 

2 2 

(s,,...,s  vary  with  C).  We  assume  b.  -i  d.  > 0,  hence  p.  , ^ p.  (i  = 1 , . . . , n) . 

in  11  1-1  1 

Much  of  the  paper  deals  with  angle-cons trained  interpolants,  in  which  case  the 


angles 


6(0)  = a,  0(s)  = B 


are  prescribed.  If  0(0)  and/or  e(s)  is  not  prescribed  the  corresponding  terminal 
of  C is  said  to  be  free,  and  the  corresponding  natural  end  conditions  for  an 
extremal  interpolant  turn  out  to  be 

(2.3)  0' (0)  = 0,  0*  (s)  = 0 . 

The  functional  which  is  made  stationary  by  an  extremal  interpolant  E is  the 
potential  energy  (or  curvature  functional) 

s 2 

(2.4)  / le' (s)l  ds  . 

0 

The  comparison  functions  are  ta)cen  from  the  Sobolev  space  2 ~ ''*1  2^*^'^' 

functions  0 : [0,SI  IR  , which  are  absolutely  continuous  and  have  derivatives 

S i 

*^2  2 2 

0'  in  L 10, S)  with  norm  {/  (0  0'  )}  . S is  a prescribed  positive  number 

^ 0 

large  enough  so  that  the  functions  in  W satisfying  conditions  (2.1)  and  (2.2) 

1 f 2 

(if  imposed)  form  a subset  with  nonempty  interior.  In  this  paper  we  do  not  deal 
with  the  existence  of  extremal  interpolants,  but  we  start  with  a )tnown  extremal  E 
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and  investigate  whether  it  is  stable  or  not.  In  this  case,  we  may  take  S = s + 6 


0 


^2  2 
/ (cosO^n  + 2sineQf,)  + ' 

^i-1 


(2.7c) 


2 2 
/ (sine^n  - 2cos0q5)  + 2d^Tj^  = 0 


Vl 


n = 1,  . . . , n . 

The  value  of  the  potential  energy  for  the  curve  is 


U(e)  = / 
0 


fJ^(t) 


t2 


s' (t) 
E 


;(t)dt  = I (1  + / 0^  . 

i=l  s.  , 


Set 

(2.9) 


= / 
s 


i-1 


I 

i=l 


Expand  (2.8)  in  powers  of  E,  using  (2.5): 


U(E) 


(2.10) 


U„  + E 


I 


i=l 


T . U. 
1 1 


+ 


2 

e 


0'C' 


n 1 n s 

2 y T . / 0'n'  + y T.u.  + / n*  + 0(e  ) . 

1=1  s.  , 1=1  0 

1-1 


since  is  a stationary  value  of  the  potential  energy,  we  must  have,  using  (2.7a), 

s . 

n 1 u , 

I I ^ cose^jh)  = 0 

i=l  s.  , '’i 

1-1 

and  this  must  be  true  for  every  n for  which  (2.7b)  holds  and  for  which 
(2.7d)  n(0)  = 0 and/or  n(s)  = 0 

if  is  angle-constrained.  From  this  one  Infers,  by  the  usual  arguments  of  the 

calculus  of  variations  (carried  out  in  detail  in  (2))  that  6^  is  continuous,  0^ 
is  continuous  between  consecutive  interpolation  nodes,  and  there  exist  constants 
f R such  that 
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(2.12) 


2ej(s) 


A . d , sin9  (s) 
11  0 


~ + A.b. 
d . 11 


cos0^ (s) 


0, 


Moreover,  conditions 
Integration  of 

0 


(2.3) 

(2.12) 


must  hold  for 
gives 

+ A . d. cos0„ (s) 
11  0 


s^,  i = 1, . . . ,n  . 

0Q  if  the  terminals  are  free. 


u. 

^ + A.b. 
d.  11 


sin0Q (s) 


0. 


I 1 J 

and  another  integration  from  s^  to  s^  shows  that  5^  = 0.  Thus, 

(2.14)  eA^(s)  + A^cos0„(s)  + A^sin0„(s)  = 0, 

0 1 0 1 0 

s.  , <s<s.,  i=l,...,n 

1-1  — — 1 


where  we  have  set 


(2.14) 


A^  = A.d.,  = -A.b.  - ^ . 

1 11  1 lid. 

1 


1 2 

To  determine  the  multipliers  X.,X.  we  use  the  fact  that  0^  and  0’  are 

11  0 0 


continuous,  hence 


(2.15)  (a^^^  - A|]cose^(s.)  t (a^^^  - A2]sin0j^(s.)  = 0, 

i = 1, . . . ,n  - 1 . 

Conditions  (2.15)  together  with  the  interpolation  conditions  (2.1)  and  end  conditions 

0^(0)  = a(or  ®q(0)  “ 0).  ~ B(or  9g(s)  = 0),  are  3n  + 1 independent  conditions 

1 2 

for  the  3n  + 1 un)inowns  A^,Aj^,s^(i  = l,...,n)  and  0^(0),  which  together  with  the 
differential  equation  (2.13)  determine  the  interpolating  elastica  0^.  There  may  be 
many  solutions  of  these  equations,  as  shown  in  (2),  but  the  distinct  solutions  are 


isolated. 

The  assumption  d^  y 0(i  « l,...,n)  was  made  only 

2 2 

The  obtained  result  remains  true  as  long  as  b^  + d^  > 0 


to  avoid  case  splitting, 
for  i = 1, . . . ,n. 
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T 


3.  St^lbiIity  Criteria, 

We  now  look  at  the  quadratic  terms  in  the  expansion  (2.10)  for  the  potential 
energy: 


(3.1) 


n 1 n 

/ n'  - 2 I T.  / 6'n'  + I T u.  + 2 / 0:c’, 

0 1=1  s.  , 1=1  0 

1-1 

s . s . 

(<^ose^)n  = ^f  (sine^)n,  u.  = / 0'^  - 

1 ^ ®i-i  ®i-i 

Using  (2.7c),  (2.12),  (2.13)  and  (2.14),  we  can  eliminate  C in  (3.1): 

2  / e'i"  - 26  • 


^i-1 


(3.2) 


s . s . 

1 1 

= -2  / e^c 


-X .d.  / (sin0^)C  • 

11^  0 

^-1 
s^ 

i/  (cose^)n^ +b^T^ 
=i-i 

1 r ^ ^.2  2 2 

= --/  0-  n - T.u.  . 


U . 

•^  + A ,b. 
d.  11 
1 


s . 
1 


^-1 


A .d. 
1 1 


- I^+X.b. 


/ (cos0q)^ 

s^ 

f/  (sin0Q)n^ + d^T^ 


®i-l 


Tlius,  since  (®q5)  - 0)  = (9q?)  + 0)  and  {®g5)  (0)  = 


S - S _ _ jj  - 

/ 0o«’  =-2/  e-  n - I x,u. 


i=l 


and  (3.1)  becomes 


(3.3) 


f (n'^  - 2 j T.  / ^ 0'n-)  - I / 6;,v  . 


i=l 


We  introduce  the  subspace  VOg)  “f 
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i * 1 , . . . , n ; 


By  f.imiliar  arqumcjnts  it  follows  that,  the  continuous  form  attains  thf?  value? 

s ^ 

for  some*  € V^(O^),  / * 1.  The  Ruler  cfjuation  for  is: 


n“(r,)  f (s)  (s)  - d.Nj  cosO^^  (s)  + cosO^n^)  0”  (s) 


♦ Pj^[b.cos0^(s)  + fl.sin0^{s)l  + li^ri^(s)  = 0, 


? , < s s. , i * 1, . . . ,n?  n! 

1- 1 — — i * 


continuous  . 


Tfj#'  multipliers  p.  € IB  result  from  the  side  conditions  / (b.cosO^  + d.sinO  ) ri  * 0, 
1 *'  1 0 i 0 

’’i-i 

s 2 

and  M * * I*'  from  the  f:ondition  / n *=  1.  It  should  be  understood  that  d^  cosG^ 

in  the  two  integral  terms  of  (3.7a)  is  replaced  by  -b^^sinO^  if  d^  = 0.  (3.7a)  is 

s»jf^pJejm*nted  by  the  conditions  of  (3,4) 


/ (b^cosO^  + d^sin0^)n^  - 0,  i = l,...,n, 


( *.7c)  n(0)  a 0 or  0^(0)  = 0 and  n^(s)  » 0 or  ni(s)  = 0 

dependinq  on  whether  is  anqle-constrained  or  free.  Besides  we  have  the  conditions 

r s , s . , -1 

-1  ^ -1 

(*.7d)  ^*1  / cosG^n*  - d^^^  / cosG^n.  =0,  i » l,...,n-l 

L 

resulting  from  the  fact  that  the  coefficient  of  n*  in  (3.5)  is  discontinuous.  If 

G* (s  ) / 0 for  the  internal  nodes  s.,...,s  , then  (3.7d)  combined  with  (3.7c) 

'/  i 1 n“l 

requires: 


(3.8)  d^^  / *^*osG^n*  ■ “^1^  / “ constant  for  i » l,...,n  . 


I 


The  multipliers  can  be  eliminated  from  (3.7a).  We  integrate  (3.7a)  over 

the  interval  (s.  ,,s,)  and  obtain: 

1-1  1 

^i  ^ ■ 

^i'^i  " ”**Vi’  ■ ''*‘®i’  ■ 2 ^ ^ ®o'’» 


®i-l 


(3.9) 


i-1 


With  (3.9)  substituted  in  (3.7a),  we  obtain  the  equation 


;(s)  + n*(s)  + 6j^(n,)cos0Q(s)  + 5^ (n*)8ine^(s)  + p*n*(s)  = o, 


(3.10) 


s^  ^ £ s ^ s^,  i + l,...,n 


where  the  6.  and  4.  are  well-defined  linear  functionals,  depending  only  on  6„. 

11  o 

(3.10)  together  with  (3.7b,c,d)  is  a nonconventional  linear  boundary-value  problem 
for  n,,  P,  being  the  eigenvalue.  Introduce  the  linear  operator  R with  domain 
D(R)  of  functions  n : 10, s]  ■*  JR,  with  h'  continuous  on  (0,sl,  n"  continuous 
on  each  (s^  l'®i'  ^ satisfying  conditions  (3.7b,c,d),  defined  by. 

(Rn)  (s)  = -n”(s)  - ■jeQ^(s)n(s)  - 6j^(n)cose^(s)  - 4^(n)sine^(s) , 


(3.11) 


3, , i - 1, . . . ,n 


'i-1  “ i' 

Then  the  above  eigenvalue  problem  may  be  stated  as: 

(3.12)  Rn  - pn  . 


A simple  calculation  shows  that  if  / n^  “ 1 then 


(3.13) 


u • / nRh  " Q(0Q»n)  . • 


Therefore,  p,  = COq^O*)  is  the  smallest  eigenvalue  of  R. 

We  conclude  that  the  form  Q is  positive  definite  if  and  only  if  p,  > 0, 


or  equivalently,  all  the  eigenvalues  of  R are  positive.  We  have  obtained 
Proposition  3.2.  The  extremal  interpolant  is  stable  if  and  only  if  the 

operator  R defined  above  has  only  positive  eigenvalues. 


The  following  proposition  provides  a useful  sufficient  condition  for  instability 
of  interpolating  elastica. 

Proposition  3.  3.  Suppose  with  normal  representation  s 9^(s)  (s^  1 ^ 

is  an  angle-constrained  extremal  interpolant  for  some  configuration  fp, .P^»..-/P  }■ 

12  n 

Suppose  E is  another  extremal  interjxjlant  (angle-constrained  or  free),  with  normal 

representation  s 0(s)  (s„  <s<'s  ,,s  <s,s  ,'-s),  where  9 is  an 

0 — — n+1  0 — 1 n+1  — n 

extension  of  0 with  no  additional  linot;  thus,  0"  is  continuous  at  s,  (s  ) if 
o In 

s„  < s,  (s  ’>  s ) . Tlien  K is  also  unstable. 

0 1 n+1  n 

Prrx)f.  The  extremal  E,  which  interpolates  the  configuration  (p.,p_,...,p  ,,p  .} 

0 2 n- 1 n+1 

can  also  be  considered  as  an  extremal  E which  interpolates  the  configuration  with 

p.  (p  ) inserted  between  p„  and  r'.,  (p  , and  p ,)  if  s„  < s,  (s  , > s ).  Let 

1 n 0 2 n-1  n+1  0 1 n+1  n 

0 denote  0 in  this  identification.  Since  E^  is  unstable  there  exists,  by 

Projxisition  3.1,  ^ snch  that  P^*'*^i^nlar , 

n (s  ) = (s  ) = 0.  Let  n be  defined  as  an  extension  of  n„: 

0 1 0 n O 

n (s)  = n_ (s)  , s ^ s < s 
0 1 — — n 

(3.14) 

=0,  s„  < s < s,  and  s < s < s 

0 — 1 n — n+1 

It  is  easily  checl<ed  that  n f Q(0,n)  = follows,  again 

by  Proposition  3.1,  that  E is  unstable.  Since  E is  obtained  from  E by  the 
removal  of  constraints,  E is  unstable. 

Let  of  Proposition  3. 3 be  angle-constrained  at  one  terminal  only,  say  at  p^. 

For  this  case  we  have  the 

Corollary.  If  the  unstable  extremal  of  Proposition  3.3  is  angle-constrained  only  at 

p and  0 is  an  extension  of  9„  to  s„  < s < s . , with  0"  continuous  at  s 

0 0 — — n+1  n 

then  E is  unstable. 

The  proof  of  this  is  an  obvious  modification  of  that  for  Proposition  3.3. 

Another  useful  sufficient  condition  is  expressed  in  the  following 
Proposition  3.4.  Supixjse  E is  an  interpolating  elastica  angle-constrained  at  none, 
one,  or  both  terminals)  and  E^  is  a subarc  between  consecutive  nodes  of  E.  If  E^, 
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f 


considerud  as  a 2-point  extremal  interjwlant  which  is  anqle-constrained  at  the 

terminals  which  are  internal  nodes  of  A',  is  unstable  then  f'  is. 

Proof . Suppose  s 0(s)  (0  ^ s <■  s)  is  the  normal  representation  of  fT  and 

s ►►  0.(s)  (s.  , < s < s.)  is  the  restriction  of  0 which  represents  R.  Since  R. 

1 1-1  — — i i 1 

is  unstable  there  exists  ^i  ^ that  Q(0j.nj^)  ^ 0.  In  particular, 


n . ( s , , ) - 0 

1 1-1 

and/or 

n.  (s. ) 
1 1 

- 0 if  i i 2 

and/or 

i 

•"  n - 1.  The  extension  n^  with 

value  0 on 

10,  s . 

1- 

i ^ 2)  and  on 

(if 

i < n - 1)  is  continuous,  and 

clearly  n t 

© 

o 

> 

Q(0,n) 

■ 1 0- 

Hence 

R 

is  unstable. 
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4.  Extremals  close  to  stciblc  ones. 


If  is  an  extremal  interjjolant  of  some  configuration  f 

s 0^(s),  0 < s < s,  is  its  normal  representation,  and  0 = s <'s  *s 

0 — — 0 I n 

are  its  interixjlation  nodes,  then  t = 0^(st)  = O^tt),  ^ ^ t ^ 1,  is  the  normal 

representation  of  an  extremal  interpolant  for  the  configuration 

{s  },  with  interpolation  nodes  0 = t_  < t < ...  < t =1, 

0 1 n 0 1 n 

t.  = sis.  If  the  terminals  of  are  free  or  angle-constrained,  so  are  those  of 

1 i ' 0 

F^,  Clearly,  F^  is  stable  if  and  only  if  F^  is.  In  the  following  we  will  often 
use  the  standardized  normal  representation  of  elastica. 

i^et  denote  the  metric  space  of  functions  &.*  (0,1]  -*■  JR,  for  which  there 

are  real  numbers  u » a(0),  B - B^O)  such  that  the  equations 

^'^(t)  = usin0(tj  - BcosO(t) 

(4.1) 

0*  (t)  * acos0(t)  + Bsin0(t),  0 ^ b £ 1, 
hold  with  the  distance  functional  ~ < t < l^^l^^^  * 

proof  of  the  proposition  below  we  will  use  the  following 


U.'nvna . 

The  functionals  fi,6  from  F^ 

to  » 

are 

continuous. 

I’roof . 

First, 

a(O),B(0)  are  uniquely 

defined 

, for 

if  (4.1)  and  also 

UjSjnO  - 

BjCosO,  0"  » Oj^cosS  + 

B^sin0 

hold, 

then 

0 - (fi  - u^,sin0(t)  - (3  - Bj^)cosO(t)  = (a  - aj^)cosG(t)  + (6  - Bj^)sin0(t),  hence 
ft  ■ and  3 * 6^.  Clearly,  (4.1)  is  equivalent  to 

1 

(4.2)  e(t)  - (1  - t)e(O)  - tO(l)  ■=  / g(t,T)  {acose(T)  + B8in6(T))dT, 

0 

whe  re 

g(t,T)  = (t  - 1)T,  0 ^ T £ t 

= (T-l)t, 

The  uniqueness  of  a, 6 In  (4.2)  Implies  that  the  continuous  functions 
1 1 

X ” / g ( • , T)cos0 (T)dT,  y * / g( • ,T)sin0 (T)dT  are  linearly  Independent. 

0 0 


I 

( 
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2 2 2 

Therefore,  the  Gramian  / x / y - (/  xy)  is  ^ 0.  Thus,  if  (4.2)  is  dot-multiplied  < 

by  X and  y respectively,  two  independent  linear  scalar  equations  for  n,6  are 

obtained,  whose  solution  demonstrates  the  assertion  of  the  Lenma. 

We  also  observe  that  the  functionals  a(9),B(9)  are  uniquely  determined  by  the 

restriction  of  9 to  any  subinterval  of  [0,11. 

Now  let  F denote  the  class  of  interpolating  elastica  E with  normal 

n I 

representation  t B(t),  0 ^ t ^ 1,  all  satisfying  the  same  type  of  end  conditions 

(free  or  eingle  constraints)  and  having  (n  + 1)  interpolation  nodes 

0=t  <t,  <...<t  =1,  where  t,  = t. (9)  for  i*l,...,n-l.  We  also  assume 

0 1 n i 1 

that  no  two  consecutive  interpolation  points  of  E coincide.  F^  is  made  into  a ‘ 

metric  space  by  use  of  the  distance  functional 

(4.4)  d(9,,9,)  = max  |t.(9,)  - t. (9  ) ( + max  |9  (t)  - 9 (t) | . 

^ ^ i=l,...,n-l  ^ ^ ^ 0<t<l  ^ ^ ' 

We  now  prove 

Proposition  4.1.  Suppose  E^  is  an  interpolating  elastica  with  normal  representation 

f , which  is  stable.  Then  there  exists  6 > 0 such  that  every  interpolating 
0 n 

elastica  E with  normal  representation  ® ® which  d(6g,9)  < 6 is  also  stable. 

Proof.  For  every  9 e F we  have  by  (2.12) 

^'^(t)  - a^(9)sine(t)  - (9)cos9  (t) , 

(4.5)  9"(t)  - o^(9)cos9(t)  + 8j^(9)sin9(t), 

t^_^(9)  ±t  < t^(9),  i “ l,...,n  . 

We  first  take  6^  > 0 so  that  d(0,0Q)  < 6^  Inplles 

(4.6)  ^ ♦ for  i - 1 n . 

It  then  follows  from  the  above  Lemma  that  we  can  find  so  that  for 

d(9,0Q)  < 6^: 

|o^(9)  - < 1.  IB^O)  - ^ ^ 


! 

■*'i( 


hence  by  (4. 5) 


max  I 0" (t) I ^ M 
O^t^l 

for  some  M.  Thus,  the  family  {0*  : 0 e /■’  ,d(O,0  ) < 6^}  is  equicontinuous : 

T\  0 2 

(4.7)  |0'  (f  ) - 0'  (t")  1 Mlf  - t"  I . 

Suppose  £ > 0 is  given.  Using  the  Lemma  again  and  Equations  (4.5),  we  can  choose 
5^>0,  that  d (0 , 6^ ) < 6 ^ implies  |u^(0)  - ^ 

is  so  small  for  i - l,...,n  that  (4.5)  yields 

(4.8)  |0' (t)  - 0^(t)  I < for  te  1 1 . (6^) , (0^)  J n ( t , (0 ) , t . (6 ) I , 

i = 1 , . . . , n . 

Let  the  overlapping  of  the  two  intervals  in  (4.8)  occur  so  that 

t.  ,(9  ) ' t.  ,(0)  < t.(0„)  < t.(0).  If  5^  > 0,  6 < 6 is  such  that 

1-1  0 — 1-1  lO-i  4 4—3 

Mit,  ,(0)  - t.  , (e„)  I < C/3  for  d(0,0„)  < fi  ^ then  by  (4.7)  and  (4.8),  for 
' 1-1  1-1  0 04 

t.  , (0„)  < t < t.  , (9)  : 

1-1  0 - — 1-1 

|0'(t)  - 9'(t)|  < Ifl'(t)  - 0'(ti.i(e„))|  t l0'(t..^(9Q))  - e'(t..^(0o))l 


7 |9'(tj.j(0o))  - 9Mt)|  ' I + I 7 ^ 


and  the  same  result  is  obtained  for  ^ . (0^)  ± t^(0).  Altogether  one  finds  that 


for  d(0,0^)  6 : 

0 4 


|9' (t)  - 0^(t) I < c,  0 < t < 1 . 


For  0 € f’  , n c W „10,1|  define  (compare  (3.5)): 
n 1,2 


1 n *^1 

(4.10)  Q(0,n)  * / (n'^  - + 2 I f (cose)n//  sin9  / en', 

0 i=l  /t..^  t..^ 


where  t.  stands  for  t. (0) . In  (4.10)  it  is  assumed  that 
1 1 


I sin9  / 0;  if  / sin9  = 0 then  the  term  in  brackets  is  to  be  replaced  by 
t.  , t.  , 


/ (sin0)n//  cos6  . It  follows  from  (4.9)  that  one  can  find,  for 


bounded  set  B < W ^(0,1),  6_  > 0,  6^  < 6 . such  that 
1,2  5 5 “ 4 

(4.11)  |Q(0,n)  - QCOg.n) 1 < 2c 

for  all  n f B and  Q € F , d(0,0  ) < 6 . 

n 0 5 

For  Q f F we  also  define  the  subspace  V (0)  of  W,  .(0,11  (see  (3.4)): 
n 0 1#  2 

t.  t. 

1 i 

v^{6)  = (n  f / dTn(t)cos(0(t)  - 0(T)l  = 0, 

’^i-i  *"1-1 

(4.12)  i = l,...,n;  and  n(0)  ■=  0 and/or  n(l)  =0  if  the 

elements  of  F are  anqle-constrained  at  the 
n 

corresponding  terminal). 

For  the  given  bounded  set  B c W |0,11  (B  is  then  totally  bounded  in  L |0,1)) 

1,2  2 

one  can  choose  5 > 0,  6 < 6 , so  that  the  L_  - distance  of  the  sets 

D 6 “ o 2 

V„(0)  n B,  V„(0„)  B is  arbitrary  small  if  d(e,0„)  6,.  Prom  this,  together  with 

0 0 0 O D 

(4.11),  one  concludes  that  6 > 0 can  be  found  such  that  d(0,0^)  < S implies 


inf  CISq.H) 


inf  Q(0,n)  < inf  Q(0Q,n) 


r^2 f^2. 


ntv  (0  ),  / n -1  n«v  (0),  / n =1  nfv^(0j^),  n =i, 

0 0 


hence  that,  by  Proposition  3.1,  £’  is  stable. 

Proposition  4,2.  If  0jj  in  Proposition  3.1  is  strongly  unsteible  (i.e. 

inf  Q(6  ,n)  < 0),  then  there  is  6 > 0 such  that  the  elastica  E for 

neV(0Q),  /n^-1  ° 

which  d(eQ,0)  < j are  also  unstable. 

We  apply  Proposition  4.1  to  extremals  which  interpolate  configurations  close 


to  the  ray  configuration.  Suppose  E^  is  the  extremal  interpolant  with  normal 

representation  0p(t)  ”0,  0 £ t £ 1,  which  interpolates  the  ray  configuration 

where  p°  - (t°,0),  0 - t°  < t°  < . . . < t°  - 1,  and  has  free 
Oln  11  ui  n 
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terminals. 


It  was  proved  in  [2,  Theorem  7.41  that,  given  e > 0,  there  exists 


6 > 0 such  that  for  every  configuration  (p  ,p  , ...,p  } with  Ip.  - p°|  < 6 

O 1 n 11 

is  a unique  extremal  interpolant  with  free  ends  and  normal  representation 

1 2 

for  which  ®q  stable.  In  fact,  Q(6.,n)  = / n'  and 

0 

*^i  1 

V^(6^)  = (n  : / n = 0,  i = l,...,n}.  In  particular,  we  must  have  / n = 0 


there 

6 

E 

for 


n c follows  that  Q(9g,P)  ^ . Thus  we  have  obtained 

Proposition  4.3.  For  every  configuration  sufficiently  close  to  the  ray  configuration  ‘ 

there  exists  a unique  stable  extremal  interpolant  with  free  terminals  that  is  close 
to  the  trivial  interpolant. 

Of  course,  this  proposition  holds,  a fortiori,  for  extremal  interpolants  with 
angle  constraints. 
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5.  Instability  of  the  2-point  interpolants  E ,E* . 

n n 

If  the  configuration  to  be  interpolated  consists  of  two  points  Pq'Pj^  then  the 
elastica  E^  has  normal  representation  0^  e C^tO.l),  satisfying  the  equations 
(see  (2.12),  (2.13)): 

- X^sin0g  + l^cos0^  “0,  ®Q  “ i^cos0Q  - A^sinS^j  • 0 . 

In  the  sequel  we  will  arrange  it  so  that  ®q  “ 0 when  ®q  ” then  these 


equations  become 


“ A„sin9_,  6"  » X cos6. 
2 0 0 0 0 0 0 


for  some  X^^  e » . if  Pq  “ (0,0),  p^  = (0,d),  d > 0,  then  the  interpolation 
conditions  are 

1 1 

f cos0  = 0,  / sine  “ d . 

0 ” 0 ° 

If  0Q(O)=a,  0^(1)  -=6,  0£a£it,  0£6;J^ii,  then  by  (5.1): 

6 0 

(5.2)  ' 1/  sin’^'^^u  du] , (2XQ)^^^d“  |/  sin^'^^u  du|  . 

a a 

However,  these  formulas  for  X^  and  d are  correct  only  if  0 0 < t < 1 

(i.e.  E^  has  no  internal  inflection  point);  otherwise  they  must  be  modified,  as  will 
be  done  below. 

The  quadratic  form  (3.5)  becomes 

1 , 2 ^2 

(5.3)  eO./h)  - / (n*  - X sin0  n ) - (2X  /d)(/  cose  n) 


and  it  is  to  be  minimized  on  the  space  (see  (3,4)); 


“ tn  * W ,(0,1I  : / nsln0  - 0; 

0 0 1,2  ' 0 


n(0)  - 0 and/or  n(l)  ”0  if  E^  is  angle-constrained)  . 


a.  We  first  Investigate  the  stability  of  the  extremal  i with  free 

terminals  which  has  (n  - 1)  internal  and  2 terminal  inflection  points.  E consists 

n 
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of  n arcs,  comiruent  to  E whicli  is  the  basic  nontrivial  2-rxjint  extremal 

interrolant  (see  [2,  Sec.  5]).  Tf  0 is  the  normal  representation  of  E and  we 

n n 

choose  0 (0)  = 0 then  0 (t)  varies  from  0 to  ti  to  0 to  . . . to  —[I  - (-•l)*^l'n 
n n 2 

as  t varies  from  0 to  l/n  to  2/n  to  ...  to  n/n.  The  points  k/n(k  = l,...,n  - 1) 

are  the  internal  inflection  points.  The  total  variation  of  0 is  Va{G  ) = nti  . 

n n 

We  have 
(5.5) 


'^)'^(t)  = X sinO  (t) 
2 ri  n n 


0 (l/n  + t)  * TI  - 0 (t),  e (2/n  + t)  0 (t) 
n n n n 


0 (0)  - 0 
n 


Formulas  (5.2)  are  now  replaced  by 


(5.6)  (2X  ) " n / sin  u du,  (2X  ) d » n / sin  u du  . 

n 1 n ' 

0 0 


We  choose 


ri  = 6 - d^J  9 sin9 


1 1 
Then  f nsinO  = 0,  hence  n e V.,(G  ).  Since  / cosO  = 0,  we  have 

b " ” " 0 " 


\ ncosO  * f 9 cosO  , and  (5.3)  becomes 
6 0 ” " 

1 1 ^ 1 2 
(5.H)  0(0  ,r|)  = 2>  d - X / S'^sinO  ♦ (X  /d)(/  9 sinO  ) - (2X  /d)  (/  0 cos0  ) 

^n  n njl^nn  n^nn  n^nn 


We  use 


1 2*2  * ^2 
(/  0 SinO  ) ‘ / 0 sinO  / sinB  ® d / 0 sinO 

o''  ""o'’  " 0 " 0 ” 


and  find 

(5.9) 


p(0  ,n)  ^ (2A  /d) (d^  - (f  0 cosO  )^)  . 

* n — n ' n n 


1 


To  evaluate  the  integral  term  in  (5.9)  we  first  assume  n even.  Then 


1 1/n  1/n 

/ 6 cos6  = ■?  / 9 cosO  - / (IT  - 6 )cos6 

0 " " ^ 0 " " 0 


= -(n/2A  )2  / 0'^  = -2d  . 


We  find  the  same  result  for  n odd.  (5.9),  (5.10)  show  Q(6  ,ri)  0.  Thus,  we 

n 

have  proved  that  E is  unstable.  This  was  also  proved  in  [21 , but  by  a different 
n 


b.  We  now  show  that  the  above  extremal  is,  for  n ^ 2,  also  unstable  if 

anqle-constrained  at  both  ends.  Let  this  extremal  be  denoted  as  E* . If  0*  is 

n n 

1 2 

its  normal  representation  then  0*  minimizes  / 6'  among  the  functions  that 

^ 0 
1 1 

satisfy  the  interpolation  conditions  / cos0  =0,  / sin0  = d and  the  end 

0 0 

conditions  0(0)  = 0,  6(1)  = -^(l  - (-D^Iti.  E*  coincides  with  E of  paragraph  a., 

2 n n 

1 

hence  6*  = 0 . n f V(9*)  now  requires  n(0)  = n(l)  =0  besides  / nsin0  = 0. 
n n n 0 ^ 


We  choose 


ri^(t)  = 9'  (t)  , 0 < t < 2/n 

• n — — 


0,  2/n  < t < 1 . 


Then,  clearly,  n f V„(9*),  and  also  / n„cos0  = 0.  Thus  (5.3)  becomes 
* 0 n 0 * 


2/n 

Q(0*,n.)  = / (n;  - A^sine^n.) 


2/n  _ „ 9 9 2 2 

= / (A'^cos'^e  - 2A‘‘sin  e ) * A^12/n  - 3 / sin  6 ). 
^nnnnn  ^ n 


But,  using  (5.5),  (5.6)  and  integration  by  parts,  we  find 


2/n  t _,,j 

(5.13)  / sin  0 = (2//2A“)  / sin  ' u du  = (2/3/2A7)  / sin  ' u du  = 2/3n  . 

0 " " 0 " 0 

Thus,  Q(h.)  * 0,  and  this  proves  instability  of  the  extremal  £*,  n ^ 2. 

In  the  next  section  it  will  be  proved  that  is  also  unstable. 


i 

I 

< 
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6.  Two-point  angle-constrained  interpolants  with  no  inflection  point. 

In  this  section  we  prove  that  2-point  angle-constrained  interpolants  are  stable 
if  they  have  no  inflection  point,  and  are  unstable  if  they  have  at  least  2 inflection 
ixjints. 

Proposition  6.1.  A 2-point  angle-constrained  extremal  interpolant  E with  no 
inflection  point  is  stable. 

Proof ■ If  E has  no  inflection  [joint  then  E is  a proper  subarc  of  the  basic 
2-ptjint  extremal  E^  (see  Sec.  5).  Clear./  E is  contained  in  another  proper 
subarc  E of  E^  which  has  an  axis  of  symmetry.  By  Proposition  3.3  it  suffices 
to  prove  that  the  angle-constrained  extremal  E is  stable.  Let  t 9(t)  (0  t < 1) 
be  the  normal  representation  of  E and  (0,0),  (0,d),  (d  > 0)  the  coordinates  of 
the  terminals,  with  0=0  along  the  positive  x axis.  Then  we  have  the  following 
equations  for  0: 

^'^(t)  = lsin0(t),  0 ^ t ^ 1 


(6.  1) 


0(t)  = It  - 0(1  - t) 

0 (0)  = a,  0 < a < it/2 

1 V -1/7 

d = / sin0  = 2(2A)  ' f sin  ' u du 


= 2 / sin'^^^u  du  . 


It  follows  from  Proposition  4.1  that  E is  steible  for  all  a sufficiently  close  to 
ii/2.  Hence,  if  E is  unstable  for  some  a > 0,  there  exists  a smallest 
a = a^,  0 < < tr/2,  for  which  E = E^  (correspondingly,  unstable. 


It  then  follows,  by  Prcoosltion  4.2,  that  inf 


n*Vo(0o).  /n  =1^'  0 


2 ,Q(0_«n)  ” 0,  hence  there 


exists  e V^(0^),  n^  ^ 0,  such  that  Q(0p,nQ)  = 0.  We  will  show  that  this  is 
not  the  case. 
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By  (3.4)  and  (3.5)  we  have 


(6.2) 


2*°'^'  ’ ^ nsinC^  . 0) 

1 - 1 . 
Q<0  ,n)  » / (n’  - A sln8  n ) - (2A  /d)  (/  cose  n)" 


^'’^nev(0  ) ° implies  (see  Proposition  3.2  and 

Equations  (3.7a,b)  that  Pp  satisfies  the  following  system  for  some  Pq  * ® • 

n'Mt)  + A sine  (t)n-(t)  + o cos9  (t)  + p sin6  (t)  » 0 
\J  uuu  uu  uu 


(6.3) 


n (0)  = ^ '’o  ° 

0 


% = ‘^Ap/d)  / n^coseg  . 

The  equation  n"  + A^sine^n  “ 0 has  the  general  solution 


(6.4) 


■ =0®0  ^ ^®0^0'  ^0*’^’  ■ I (l/8ineg(T))dT 


By  using  the  method  of  variation  of  parameters  one  finds  for  the  general  solution  of 
the  differential  equation  in  (6.3): 

(6.5)  ng(t)  . -(o^/2A^)t9'(t)  - (p„/A^)  ^ c^e'(t)  ^ c^e-(t)y^(t)  . 


0^(0)  » n^d)  ” 0 give,  since  9g(0)  » ®o(i)  =“  •'o 


(6.6) 


'=0  “ *=1  ■ • 


By  the  use  of  integration  by  parts  one  finds 


(6.7) 


/ PjjCosejj  “ - (Op/2Ajj)  (sinOp  - d)  + Cj^(^(l)sinaQ  - 1) 


and  since  this  must  equal  (daQ/2A^)  by  (6.3),  one  obtains 


(6.8) 


■ (Oo/2Ao)sinao  + Cj^(Yq(1)  sina^  - 1)  - 0 


(6.8)  together  with  (6.6)  gives  "0,  Th'iSf  we  are  left  with 


(6.9) 


Pq  “ ‘Po''^o''o>  % ■•'o’ 
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The  final  condition  / n^sinO^  = 0 yields 

(p  /A  < ) (2cosa„  - K„d)  = 0 . 

O U 0 (J  0 

since  * 0 implies  * 0,  we  must  have 

0 ■ • * cosOq  - * 

By  (6.1)  wo  have  = 0^(0)  = (2AgSinaQ)  d/2  = (21^)"^''^  / ^ sin^^'^u  du, 

“o 

hence  G(0)  = 1,  G(ii/2)  = 0,  G’(a)  = -<l/2)sin  ^^^acosa  J sin^^^u  < 0 for 

a 

0 < a < 11/2.  Therefore  no  '^q'Pq  satisfying  (6.3)  exist,  and  the  proof  of 
Proposition  6.1  is  complete. 

We  prove  next: 

Propos i tion  6.2.  A 2-point  extremal  interpolant  E (angle-constrained  or  free)  with 
2 or  more  inflection  points  is  unstable. 

Proof.  If  E has  at  least  2 inflection  points  then  E contains  the  basic  2-point 
extremal  E^  (see  Sec.  5).  By  Proposition  3.3  it  suffices  to  prove  that  which 

is  E^  with  angle-constraint,  is  unstable.  We  do  this  by  exhibiting  f '^0**^1*' 

/ 0,  for  which  Q(9j^,nj^)  = 0.  As  in  the  preceding  proof,  this  will  lie  the  case 
if  for  some  ® = 

n'^  + Aj^nj^sin0j^  + Oj^cosSj^  + pj^sinBj^  • 0 
1 

(6.12)  n (0)  - n (1)  = 0,  / n.sine  - o,  n,  o 

1 1 Q i 1 1 

1 

o * (2X  /d)  / n.cose  . 

i 1 Q i 1 

This  system  is  satisfied  by 

p^  - 0,  nj(t)  - (1  - 2t)0^(t)  . 
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7.  Two-[)oint  antjle-constrained  interi>olants  with  one  inflection  point. 

If  tho  2-ix)int  anqle-constrained  extremal  P’  contains  one  inflection  point 
(either  at  one  end  or  internally)  then  the  problem  of  stcibility  is  more  complex.  If 
one  proceeds  from  the  inflection  point  0 along  t'  in  one  or  the  other  direction 
to  a terminal  one  traverses  a proper  subarc  of  the  basic  extremal  h’^  (see  Sec.  S) . 
■niere  is  a iK>int  on  close  to  the  far  terminal,  - its  precise  location  is  given 

below  - on  which  the  stability  of  f"'  depends.  We  call  this  |K)int  a stability  focus. 

E'  may  contain  the  right,  the  left  or  neitlier  stability  focus.  We  prove 
Proix>sitiQn  7.1.  A 2-ix)int  angle- cons  trained  extremal  inter[x>lant  K with  one 
inflection  x^int  is  stable  if  E contains  no  stability  focus. 

Proof.  E contains  neither  stability  focus  as  one  proceeds  from  the  inflection  point 
to  one  or  the  other  terminal,  hence  is  a subarc  of  another  extremal  E,  which  is 
symmetric  with  respect  to  the  inflection  {joint  and  also  contains  no  stability  focus. 

By  Proposition  3. 3 it  suffices  to  prove  that  the  angle-constrained  extremal  E is 
stable.  Let  t ^ 0(t) (0  £ t ^ 1)  be  the  normal  representation  of  E and 
(-b/2,  -d/2),  (b/2,  d/2),  (b  > 0,  d > 0)  the  coordinates  of  the  terminals,  with  0*0 
along  the  {positive  x axis.  We  then  have: 

|6'^(t)  = Xsin6(t),  0 < t ^ 1 
0(t)  = 6(1  - t) 

0(0)  = a,  0 < a < II,  0(1/2)  = 8’  (1/2)  = 0 


(7.1) 


b = / cos6  ' 2(21)  / cosu-sin  du 

0 0 


1 

d - / sinO  - 2(21) 
0 


a 


g du 


(21)*^^^  - 2 I sin“^^^u  du  . 

0 

It  follows  from  Proposition  4.2  that  fc'  is  stable  for  all  a sufficiently  small. 
Further  if  a ” ii,  H contains  2 inflection  points,  hence  is  unstable.  Thus  there 
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is  a smallest  a = a^,  o < < n,  for  which  E ^ F,^  (correspondingly, 

is  un.'itable.  As  one  proceeds  along  this  from  the  inflection  point  to  one  of 

the  terminals  one  reaches  the  (left  or  right)  stability  focus,  mentioned  in  the 
statement  of  the  proposition. 

By  Proposition  4.1,  we  are  left  to  find  and  0^,  so  that 

'"Vv„(0.),  /n"=l2'®‘'"^  ' 

where  § = 9^  satisfies  (7.1),  with  a replaced  by  a^.  By  (3.4)  and  (3.5)  we  have 

1 

V (0.)  = (n  e W _[0,1)  ! / n(b  COS0.  + d sin0.)  = 0} 

0 ’ 1,2  0 

(7.  3) 

1 , 1 
Q(0*,n)  = / (n'  - X^sin0,n  ) - (2Xyd)(/  ncos0,)  . 

0 0 

The  infimum  0 of  Q(0,,n)  is  attained  for  n = e if  (see  Equations  (3.7a,b)) 

Hj,  satisfies  the  following  system  for  some  P„  f 3R  : 

0*  o*cos0|,  + p|,(b  cosQ^  + d sin0^)  = 0, 

1 

(7.4)  n„(0)  = n^d)  =0,  / (b  cos0^  + d sin0^)  “0,  n,  o, 

0 

1 

d*  = (2X^/d)  / n*cos9,  . 

0 

Using  the  general  solution 

n(t)  = - [(o.  + p,b)/2X,lt0;(t)  - P»d/X,  + Co9;(t)  + CjY.(t) 


0;(t)  / (l/sin0,(T))dT  for  0 < t < 1/2 
0 ~ 

2 for  t - 1/2 


i 0;(t)  / (l/sin0^(T))dT  for  1/2  < t < 1 
t 

of  the  differential  equation  in  (7.4),  one  finds  after  lengthy  calculations, 
(7.6)  n,(t)  - (1  - 2t)0;(t)  - 0;(O),  p*  - x^e;(o)/d  . 


! 


T" 


r 


( 


Usimj  integration  by  parts  and  the  relations,  following  from  (7.1): 

2 

(7.7)  2sina^  = b = "2k^/A^,  where  = 0J,(O) 

one  obtains 

1 

(7.8)  = (2Aj,/d)  / n„CO30,,  = - A„K,,b/d  . 

0 

Then  one  verifies  readily  that  (7.6)  solves  the  differential  equation  in  (7.4);  also 

1 

1^(0)  = n,(l)  = 0 and  / n„sin0|j  = 2cosa|^  - 2b  - ic^d,  hence 
0 

(7.9)  j n„(b  cos0^  + d sin0,,)  = -2<^/X^  + 2d  cosa^  - ir„d^  . 

0 

TTius,  all  the  conditions  of  (7.4)  are  satisfied  if  the  quantity  (7.9)  is  0,  or 
using  (7.1)  and  (7.7)  and  the  abbreviation 

“ 1/2 

S(a)  = / sin  u du,  0 ^ a ^ 

0 

(7.10)  F (a^)  ; = sin^’^^Oj^S^  (o^)  + cosa^S(a^)  + 2sin^'^^a^=  0 . 

'1,,  is  the  unique  root  between  it/2  and  it  of  (7.10).  Since  F(ii/2)  > 0 and 
F(ti)  < 0,  there  is  a root  between  it/2  and  it,  and  since  F' (a)  < 0,  the  root  is 
unique  (a  rough  estimate  shows  a * 171°). 

We  have  shown  that  §,  given  by  (7.1),  with  a < a,,  is  stable  and  this  completes 
the  proof  of  Proposition  7.1. 

The  result  in  Proposition  7.1  is  sharp  because  we  have 
Proposition  7.2.  A 2-point  angle-constrained  extremal  interpolant  E with  one 
inflection  point  is  unstable  if  E contains  the  two  stability  foci. 

Proof.  In  the  proof  of  Proposition  7.1  it  was  seen  that  the  extremal  E^  whose 
terminals  are  the  stability  foci  is  unstable.  By  Proposition  3.3  E,  which  contains 
E^,  is  unstable. 


f 


I 
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There  remains  the  case  where  the  2-[x>int  angle-constrained  interjxjlant  f’’ 
contains  one  inflection  point  and  one  stability  focus.  We  may  assume  that  the  normal 
representation  0 of  R is  a solution  of 


(7.11) 


|•[e■  (t)  1^  = AsinO  (t)  , ^ t ^ 1 

0(0)  = a,  0(1)  = B 


for  some  At  IR  , where 


(7.12) 


0<a<a^£B<7i 

O(t^)  = O'(t^)  = 0 for  a unique  t^  . 


The  numbers 


and 


are  detentiined  from  the  relations 


(7.13) 


a ^ a 

/jr  = / sin  du  + ■^  sin  ^^^u  du,  du  . 


It  is  seen  that  t^  £ 1/2  and 


(7.14)  0(tjj  - T)  = Oltg  + T),  0 T < tp  . 

We  now  show  that  for  each  a,  0 < a < a^,  there  exists  a unique  S = B^(a) 
such  that  the  extremal  R is  stcible  if  6 < fi,(oi)  and  is  unstable  if  B ^ *^*^“*' 

We  say,  the  point  on  R for  which  0 has  the  value  B,,(a)  is  conjugate  to  the 
point  for  which  6 has  the  value  a.  As  a approaches  otj,  (from  below)  B*(a) 
approaches  from  eibove,  hence  the  stability  foci  of  Proposition  7.1  are  the 

special  case  of  conjugate  points  where  6*  (a)  = a. 

We  now  prove 

Proposition  7.3.  Suppose  R is  an  angle-constrained  extremal  2-point  interpolant 
with  normal  representation  t ►►  6(t),  ^ 1.  t f.  which  contains  one  inflection 
point  and  for  which  0(0)  ” a,  0 £ a a„  (see  (7.10)),  0(1)  “ B ^ a*-  ^ is  stable 

if  and  only  if  B < B*(a)<  where  B*(o)  is  the  unique  root  between  a,  and  x of 
Equations  (7.17)  below.  As  a increases  from  0 to  a^,  B,(a)  strictly  decreases 

from  X to  a^. 
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Proof.  By  ProfAjsitions  6.2  and  7.1  i'  is  stable  if  6 < and  unstable  if  B • it  . 

liut  B,(u)  denote  the  smallest  value  of  B for  which  E is  unstable  and  let  0^ 
be  the  normal  representation  of  the  extremal  /■’,  for  which  0,(0)  * n,  0,(1)  = B,(<»). 
There  must  then  exist  n,  f '^0*****'  / ^ that 

(7.15)  ^"^eV^(0,),/n'=l'^‘'’*'''’  = »'“*''>*>  “ • 

As  in  the  proof  of  Proposition  7.1,  we  have  for  n,  the  system  (7.4).  The  general 
solution  of  the  differential  equation  in  (7.4)  is  given  by  (7.5).  One  computes,  using 
integration  by  parts, 

1 1 

/ Y,cos0,  = 1,  / Y,sin0,  = 2/9'^(l)  - 2/0J,(O) 

0 0 

(7.1b) 

1 1 
/ t COS0,  (t) -OJ,  (t)dt  = sinB  - d,  / t sin0,(t) -0' (t)dt  = b - cosB  . 

0 0 

1 

The  four  conditions  n,(0)  « n,(l)  =0,  / n,cos0,  = do,/2X, 

0 

1 

/ (b  COS0,  + d sin0,)n,  = 0 for  n,  £ the  condition  n,  ^ 0,  then  lead 

0 

to  the  equation 

1/2  2 

H(a,B)  :=  (sinasinB)  (S(a)  + S(B)) 

(7.17)  + (sin^^^acosB  + sin^'^^Bcosa)  (S (a)  + S(6)) 

+ 2 (sinasinB)  (sin^'^^a  + sin^'^^B)^  = 0 

for  8 “ B, (a) . One  finds  that  the  function  B H(a,B)  is  strictly  decreasing  for 
a,  ^ B Ti.  Also,  if  0 < a < a,, 

-1/2 

(7.18)  H(a,iT)  < 0,  H(o,a)  •=  4F(a)sln  ' a > 0, 

where  F in  the  function  of  (7.10)  and  F(a)  > 0 since  a < a,.  It  follows  that 
B,(a)  is  uniquely  defined  by  (7.17).  Then  0 =•  0,,  with  0,(0)  » a,  0,(1)  = B,  (a)  , 
is  the  normal  representation  of  an  extremal  £,,  for  which  there  exists  n,  e 
with  n,  “ 0>  such  that  (7.15)  holds.  Therefore,  ff,  is  unstable,  and  by 
Proposition  3.3,  £’  is  unstable  if  E contains  £,,  i.e.  if  6 ^ B,(a). 


The  function  a ll,(a)  is  non  increasing.  For  if  ’^2  ^ '^l'  * 

then  the  angle-constrained  unstable  extremal  with  6^(0)  = a^,  6^(1)  = (ij,(Uj), 

is  contained  in  the  extremal  E^  with  6^(0)  = a^,  which  is  stable 

since  ^ contradiction  to  Proposition  3.3.  Actually, 

is  strictly  decreasing,  for  if  B.('i  ) = for  a > a then  B,,(a)  is  constant 

tor  ^ 1 which  is  imjx>ssible  since  the  function  6„  is  analytic.  Clearly, 

B,(0)  • t and  = a,,  thus  the  proposition  is  completely  proved. 

We  proceed  to  give  a geomf;trlc  interpretation  of  conjugate  £X>ints  on  a 
simple  elastica  curve.  At  the  same  time  we  obtain  the  precise  range  of  angles  that 

an  arc  of  the  elastica,  which  contains  one  inflection  jxjint,  can  make  with  the  chord  ' 

connecting  the  endrKjints.^ 

Let  E be  the  simple  elastica  of  Proposition  7.3,  t »*  6(t)  (0  ^ t ^ 1)  its 
normal  representation,  6(0)  * a,  6(1)  * B,  with  0 ' a ^ ft  < ir , and  let  Pq»Pj^ 
the  local  vectors  of  the  terminals  of  A'.  In  the  original  Interpolation  problem  the 

length  of  the  vector  p^  - p^  and  the  angles  A,B  that  E makes  with  p^^  - p^  at 


the  endpoint  are 

prescribed. 

More  precisely,  let 

A, B denote 

the  angles  in  (-ii,it) 

from  the  vector 

Pi  - Pq  to 

the  oriented  curve  A’ 

at  Pq,Pj, 

respectively.  We 

investif^ate  the 

relationship 

between  a,B  and  A,B 

. Clearly, 

(7.19a) 

(1  - B - A - B 

If  the  inflection  [xjlnt  ir  taken  as  the  origin  of  a cartesian  coordinate  system 
xy,  with  the  positive  x-axis  along  6 ■■  0,  then  the  point  (t,e(t))  on  E has 
coordinates 

X “ / cos6(T)dT  • (2/»^2jr)  Bin*^^^6  (t)  sgn  (t  - t^) 

S 

t 

y - / 8in6(T)dt  - (l/i/2X)s (6 (t)  sgn(t  “ bp)  • 

Expressing  the  slope  of  the  vector  Pj  - p^,  we  obtain 

(7.19b)  (S(a)  + S(B))/(28ln^''^a  + 2sin*'^^Bl  - tan(a  - A)  . 

^I  wish  to  express  my  gratitude  to  Dr.  D.  D.  Pence  for  his  valuable  assistance  with 
this  investigation.  - M.G. 
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since  is  above  and  to  the  right  of  it  follows  that 

(7.19c)  0 < a - A Ti/2 

For  each  pair  (a,B),  with  atB>0,  there  is  a unique  pair  (A,B) 

with  -n  < A ^ B < n determined  by  Equations  (7.19a,b,c)  (actually  A ^ -ii/2)  . 

Let  B(A;a,B)  be  the  angle  B for  fixed  A,a,B,  and  set 

(7.20)  B*(A)  = sup  B(A;o,B)  = B(A;a  ,B  ) . 

0£a<B^Tr  * * 

It  is  readily  found  that  if  o = 0 then  A = -ii/2  and  B = and  if  B = then 

A A A 

A » -n/2  and  a,  = 0;  also  if  a,  = B.  then  a = (solution  of  (7.10))  and 

A A A A * 

A “ B*(A).  We  write 

(7.21)  A*  = B(A*;a^,aj,)  = B(A*;a,B) 

O^u^^B^n 

(A*  * 99.5°).  It  follows  that  if  A is  neither  -ti/2  nor  A*  then  the  supremum 

in  (7.20)  is  attained  in  the  interior  of  the  region  0 < a < B < ir.  Tbus,  (“.»B.) 

— — — A A 

make  B«A-a+B  a maximum  under  the  side  condition  (7.19b).  It  follows  that 

a >■  a , B > B satisfy  the  equations 
A A 

1/2  1/2 

1 + MO/3a)tS(a)  + S(B)  - 2 tan(a  - A)  (sin  ' a + sin  ' B)  1 = 0 

(7.22)  . <2  1/2 

-1  + pO/3B)(S(a)  + .5(B)  - 2 tan(a  - A)  (sin  ' a + sin  ' 6)1  = 0 . 

Elimination  of  the  multiplier  p,  and  use  of  (7.19b)  yield 

(7.23)  “‘“a'®a’  “ ° 

where  H is  the  function  (7.13).  Tbus  sup  B(A;a,B)  is  attained  for  conjugate 

values  a ,6  . This  is  also  true  in  the  excluded  cases  A = -ti/2,  A = A*  since 
A A 

(O.n),  (o,,a^)  are  conjugate  pairs.  Each  conjugate  pair  (o,B)  occurs  in  this 

characterization;  for  if  a, 6 are  used  in  (7.19b,c)  a unique  A is  obtained  for 

which  a - a,,  6 - B..  We  have  proved 
A A 

Proposition  7.4.  Suppose  A,  -x  < A < n,  is  such  that  there  exists  a simple  elastics 
t'  with  terminals  Pq/Pj*  which  contains  an  inflection  point  and  which  makes  the 
angle  A with  the  vector  - p^  at  p^.  Then  the  largest  angle  B that  E can 
make  with  p^  - p^  at  Pj^  is  obtained  If  P^/P^  are  conjugate  points  of  E. 
Conversely,  each  paii  conjugate  points  Is  characterized  in  this  way. 
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1 


Wa  proceed  to  determine  the  range  of  angles  A,li  that  a simple  elastica  with 

one  inflection  point  can  make  with  the  chord  joining  the  endpoints.  Because  of 

symmetry  it  suffices  to  determine  the  half  where  A < B,  which  we  denote  as  R 

— 

If  0 ^ A ^ A*  (see  (7.21))  then  the  interval  {A,  A < B < B* (A) } is  in  R 

■—  — A^B 

(B*  (A)  as  in  (7.20)).  If  A > A*  then  there  is  no  B ^ A such  that  (A,B)  e 

this  follows  from  the  above  discussion.  Let  us  assume  now  A < 0.  By  (7.19c),  we 

have  A ^ -ii/2:  so  fix  A,  0 < A ^ ~ii/2.  Substitute  B - A + a for  B in  (7.19b), 
which  then  defines  B as  a function  of  a.  It  is  easily  found  that  3B/3a  at  a = 0 

is  B taltes  on  its  minimum  B^(a)  for  a = 0,  hence  by  (7.19a,b) 

(7.24)  S(B^(A)  - A)  = 2 sin^'^^  (B^  (A)  - A)tan(-A)  . 

It  is  easy  to  see  that  the  interval  {A,  B^(A)  < A < B*(A)}  is  in  In  summary, 

we  have 


= {-ti/2  < A < 0;  B (A)  < B < B*  (A)  } u {0  < A < A*:  A < B < B*  (A)  } . 

■ — — — — — — 

Remar)c.  The  general  (A,B)  e U /?_  is  the  image  of  two  pairs  (a,B),  hence 

A^B  B^A 

arises  from  two  distinct  simple  elastica  E . If  angle-constrained,  no  more  than 

( A , B ) 

one  of  these  is  stable.  There  may  be  no  stable  elastica  at  all  for 

(A,B)  « Thus,  if  0 < A < A*,  B = B*  (A)  , then  B = B(A;  a, ,6.)  and 

A<B  B<A  — — A A 

there  is  a unique  E,  .,,  , whose  terminals  are  at  the  conjugate  a ,B..  By 
(A,B)  A A 

Proposition  7.3,  the  angle-constrained  E is  not  stable.  It  seems  probable 

( A,  B) 

that  this  happens  only  on  the  boundary  of 

The  last  proposition  of  this  section  deals  with  2-point  interpolants  with  angle 
constraint  at  only  one  end. 

Proposition  7.5.  A 2-point  extremal  interpolant  E which  is  angle-constrained  at 
one  terminal  and  free  at  the  other  is  stable  if  and  only  if  E contains  no  stability 
focus. 

Proof.  For  the  normal  representation  t 0(t)  (0  ^ t ^ 1)  of  E we  may  assume 


i O'^(t)  = AsinO(t),  0 ^ t ^ 1 


(7.25) 

0(0)  ==0’  (0)  = 0; 

0(1)  * 0>O  prescribed 

1 

1 

/ COS0  = b and 

/ sin0  » d prescribed 

0 0 


If  H is  suf f iclfcntly  small  then  E is  clearly  stable.  If  E is  stable  for  some 
• 0 then,  by  the  Corollary  to  Proposition  3.4,  E is  stable  for  each  fS  < 

<)n  the  other  hand,  E is  not  stable  if  B = n since  in  this  case  E is  unstable 
even  if  am) le-constrairied . It  follows  that  there  exists  B^,  0 < B,  < such  that 

E is  stable  for  B < B,,  but  unstable  for  B > B^.  By  the  same  arguments  as  in  the 
earlier  part  of  this  section  we  conclude  that  we  have 
(7.26) 


where 


I 

(7.27)  V (6)  = (n  € w (0,11  : n<l)  = O,  j r\(b  cosB  + d sinB)  = 0)  . 

0 1,2  g 

For  n,  we  have  the  conditions  (compare  (7.4)): 

n"  + An^sin6  + ocosO  + p (b  cos0  + d sin0)  « 0 

(7.23)  j 

0 = (2i/d)  / n.cosB,  n;(0)  =0,  n.  (1)  = 0,  n.  j'  o . 

0 

The  condition  hi(0)  = 0 results  from  the  fact  that  if  0 " minimizes  Q(0,n) 

then  must  satisfy  the  free  boundary  condition  nJ^(O)  = 0.  Proceeding  as  in  the 

proof  of  Proposition  7.1,  one  finds 
(7.29) 

is  a solution  provided  0 (which  enters  (7.28)  through  A, 

(2X)'^^^  * / sin  ^^^u  du)  is  a zero  of  F,  cf.  (7.10).  Thus,  6^  « a^,  the 
0 


previously  found  stability  focus. 
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H.  The  stability  tunf.tion 

In  the  two  preceding  sections  the  stat^ility  problem  was  settled  for  all  extremal 
2-j/<-jint  interj^olants.  ly.-t  t.*  now  be  an  extremal,  interpolating  a general  (n  + 1)- 
r>oints  conf iqurat iori  (p  K ^rid  free  at  the  terminals  p,  »p  . For  ease  of 

formulation  we  introduce  the 

I>e f mi •;  i ' -r. . A Mibarc  /''*  of  A’*  betweeri  two  consecutive  interior  nodes  p.  wp. 
-1  1-1  1 

(2  ^ i n - Ij  is  said  to  be  proper  if  contains  no  pair  of  conjugate  points. 

The  terrfiinal  arcs  A*  and  h*  are  proper  if  they  contain  no  stability  focus. 

1 n 

By  Propositions  i.A,  7.3  and  7.5,  A*  iJi  unstable  if  any  of  the  subarcs  A'J 
is  not  proper.  We  state  this  important  result  as 

Proj[xjs i tion  8.1.  A necessary  condition  for  stability  of  an  extremal  interpolant 

with  free  terminals  is  that  each  arc  between  consecutive  interpolation  nodes  be  proper. 

It  shoul<l  be  observed  tfiat  by  assuminq  all  arcs  are  proper  we  do  not  exclude  the 

i^resence  of  inflection  {>oints.  However  we  will  exclude,  with  little  loss  of 

generality,  inflection  j^oints  at  the  knots.  We  say  A’*  is  de comfjosabl e if  p^  for 

some  m between  1 ar»d  n - 1 is  an  inflection  jioint,  otherwise  A’*  is 

indecomijosable.  If  A*  is  decomt>osaljle  then  the  s^ifiarcs  A’  from  p^  to  p 
* a 0 m 

and  A'  from  p to  are  (free)  extremal  interpolants,  and  it  is  readily  seen 

b m n 

that  A*  is  stable  or  unstable  if  both  A'  and  A'  are  stcible  or  unstable, 

a b 

re5ipectively  (the*  case  where  one  of  the  A’  ,A'  is  stable,  the  other  unstable,  is 
* a b 

omitted) . 

For  indecomp<^jr;able  extremal  interfxilants  H*  which  satisfy  the  necessary 

condition  of  Proposition  B.l  we  find  a computable  function  U*  of  n - 1 variables 

(n  1 the  ntjjnber  of  interi>olation  nodes)  with  the  firoporty  that  A’*  is  stable  if 

and  only  if  U*  has  a local  minimum  at  the  critical  fxiint  corresfxindinq  to  h'* . 

Let  s P 0*(s)  (0  ^ s <*  s*)  be  the  normal  representation  of  A’*,  with  inter- 
— — n 

rjolation  nodes  0 * s*  s?  < ...  s*.  Then  for  uniciuely  defined 

0 1 n In 

U? » • • • we  have 
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0**'(s)  ♦ A*sin0*<s)  u*ro80*(s) 


0 


1 2 

— 0**  (s)  - Xfcos0*(s)  ' u*sin0*(B)  » 0,  s*  • s s* 

2 i 1 i-l  - - i 


/ cosO*  *'■  b,  , / sinO*  * d.,  i * 1,.. 


0*'(0)  * 0,  0**(s*)  * 0 . 


In  addition  to  (8.1)  wc  have  the  corner  conditions 

(8.2)  0**(s*  - 0)  = e*»(s*  +0),  i - l,...,n  - 1 

1 1 

The  f>otential  energy  for  E*  is 


8.3  ^ ! 0*  = 5^  2(>.ni.  + ufd.)  . 

0 i i i 11 

0 1*1 

We  choose  an  arbitrary  number  6 > 0,  set  s*  + •S  * S,  and  extend  0*  to  the 

n 

interval  fO,S|  by  settinq  9*(s)  = Q*(s*)  for  s*  <■  s < S.  Every  function  <1  in 

n n — 

this  section  is  in  the  space  ''*1  2 constant  on  some  interval 

Is  ,S1,  where  0<s  =s(0)<S. 

n n n 

As  stated  above,  we  assime  each  subarc  /5’*  ( i = 1 , . . . , n)  of  E*  is  proper  and 

also  that  0* ' (s , ) ^0  for  i=l,...,n-l.  We  set  0*(sM“cif(i  = O,l,...,n). 

1 i i 

For  every  (n  - 1) -tuple  a * (a,, ,a  ,)  sufficiently  close  to  a*  » (a*,..., a*  ) 

1 n-1  i n-1 


the  system 


0''(s)  + A^sin0(s)  - p^cos0(s)  0, 

j O'^(s)  - Xj^cos0(3)  - u^sinO(s)  = 0,  I.  ’’j 


i i 

J COS0  = bj^,  / sin0  “ <3^,  i • l,...,n 


O’  (0)  « 0,  0*  (s  ) » 0, 

n 


0(s^)  - 


1-1 n-1 


h.iK  .1  uni'i'j*'  solution  0 f W,  ,,  with  A ,ii  # , 0 » s 'a  ' h, 

I * ^ 1 j J ri 

in  fi  sijf f iciontiy  fimail  prodsaiqnod  noiffhborhcxjcl  of  0*.  This  fonf)W5;  rt’adi  Jy  from 
ttte  fact  that  each  of  the  arcis  is  projier.  (H.4)  is  system  (H.I)  with  additional 

coiuiitions  0{sj  - replacinq  the  ronditions  (^J.2).  We  let  0 denote  the 

solution  of  (8.4),  the  ^ P^f  1'^  / - - - # p^  1 intorpolant  represented  by  Tlie 

jK>tential  enerqy  for  A’  is 


u„(f:  ) = / n-''  - )’  2(A,b.  + g.d.) 

On  * n /\  1 i j 1 

0 1-1 


We  now  introduce  the  function 
(8.6)  U*(rt)  = 

and  call  it  the  stability  function  (associated  with  the  extremal  /'/*).  It  is  defined 
in  a neiqhborhocxJ  of  a*.  We  prove 

i*roi>OHi  tion  8.2.  There  is  a neighborhood  N(ni*)  C ih  ^ ^ of  a*  such  that  a*  is 
the  unique  critical  i>oint  in  N(f**)  of  the  function  ti*. 

Proof.  let  W^  2 denote  the  metric  space  of  functions  ® ^ 2 intcri)olato 

the  i>oints  p.  at  nodes  s.  »s.(0),  A = s <s,  <...<s  <s*,  with  the  metric 

i 1 0 1 n 


(R.7)  d (0  ,0  ) * max  |s.(f).)  - s.(0  )|  + |0  (0)  - 0 (0)  | + ^ 

^ ^ 1 ^ 1 ^ Iq 


ir  - 'i’l' 


We  can  choose  6*  > 0 so  that  the  following  throe  conditions  are  satisfied:  (i)  0* 


is  the  only  extremal  in 


(ii)  for  each  a in 


N(9»)  = {0  E ‘*12'  1 


(8.8ii) 


N(a*)  = {a  c K : (a  - a* I < 5*} 


system  (8.4)  has  a unique  solution  0 e N(0*)  and  each  restriction  8 i , , 

a a I I®  J I 

(i  = i,...,n)  is  proper;  (ill)  for  j » - 1 


(8.8iii) 


sgn9’(s.  - 0)  » sqn0'(s  + 0)  ” sgn0*'(sj)  . 
a j a j J 


To  prove  the  proposition  it  suffices  to  show  that  a e N(a*)  is  a critical  point  of 


U*(a)  if  and  only  if  a n*. 
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By  (B.5)  we  have 


U*(a)  = y (X.b.  + n.d  ) 
i=l  ^ ^ ^ ' 

whore  the  ij^  = ij^{a)  are  determined  from  the  interpolation  and  end 


conditions : 


3 . («  . , / ot . , A , , p . ) : * / COS0  - b . a 0 

1 1-1  111  * a i 


,,a.,A,,iJ.)  :*  / sinO  - d.  » 0,  i « l-...,n 
1 1-1  111  ' a i 


E (n  ,A  ,M  ) :*  A cosa  + p sina  » 0 . 
nnnn  n n n n 

Wo  now  seek  critical  points  of  U*  as  a function  of  a ■ (a  ) and  the 

1 n-1 

accessory  variables  A -.,,,A  , p , ..,,p  , a / a , under  the  2n  + 2 side  conditions 
1 n 1 n 0 n 

(8.10).  If  a is  a critical  point  then  there  exist  multipliers  p^,o^(i  « l,...,n) 

and  such  that 

I n 


-i-  / y 


‘^‘’k  ^ ^ "k^k  * '’k‘’k>  * “ 


,E,  + li)  E V 
11  n n) 


where  Y stands  for  each  of  the  variables  • 

Let  first  i (2  ^ i ^ n - 1)  be  such  that  E*  has  no  inflection  point.  Then 
Stjn0»'(s*  ) - sqne*'(s*)  = 1,  say,  and,  by  (8.8iii),  0'(s)  >0  for  s,  , < s < s, . 

1-1  1 (j  i-1  — — i 

Thus,  using  (8.4),  we  find 

“l 

®1 f“i_l (u)c«su  du  - b^ 


*^i  *”'i-l'“i'^i''‘i^  “ / (u)sinu  du  - d^^. 


Kj^(u)  - {2X^cosu  + 2u^slnu)^ 
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Usinq  ^ \ = U j in  (8.11),  one  obtains 


b.  - p / 


1 X ' 1 

“i-l 


3 2 f 

, cos  - o.  J K,  sin -cos  = 0 


.-3_,. 
i 


i-1 


a , a . 

/■  ^ -3  f ^ -3  , 2 „ 

d.  - p.  / K.  cos-sin  -0.1  < , sin  = 0 

1 1 1 1 ' 1 
a.  , , 

1-1  1-1 


since  b.  = / cosO  = f k , ^cos  = / k,^(2X  cos  + 2p,sin)cos  : 

i'  (I'l  x-i  X 


-3, 


1-1 


i-1 


i-1 


i _3  2 ^ -3 

(2A.  - p.)  / K.  cos  + (2p,  - o.)  / K . cos -sin  = 0 
11^  1 ii  1 

a,  . a.  , 

1-1  1-1 


(8.14) 


(2A.  - p.)  / K.^cos-sin  + 
11'  1 

"i-l 


(2p  . - o , ) / k:  . 
1 1 ' 1 


-3.2 
sin  = 0 


By  the  Scliwarz  inequality 


-3 

K , cos- Sin 


2 f ^ -3  2 r ^ -3.2 

<1  K . COS  I K . Sin 

•'  1 •'  1 

"i-1  "i-1 


(equality  cannot  hold),  hence  (8.14)  gives 


(8.15) 


= 2A.,  o^  = 2p.  . 


If  i « 1 then  0**(s*)  / 0 (since  is  proper),  say  6**  (s*)  > 0,  and 

also  0*(s)  > 0 for  0 < s < s, , hence  (8.12)  holds  for  i **  1 (the  integrals 
a — 1 

involved  are  improper).  To  avoid  the  divergent  integrals  in  (8.14),  we  set 
(B.16i)  - X^G^ 
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where 


I ic^^(Aj^cob  + Uj^sin)  - - n^dj 


(8. 1611) 


1 ,,2  ^ 2,-1 

2 <\  ^ *'l> 


/ - A^bj  - Pjdj 


2 2-1 
^ ^ ‘'i> 


/ Kj^^C-A^sln  + Pj^cos)  + A^dj^  - u^b^ 
•^0 


' * ►'i>  ^*''i*“l^  * ^l^^l  ■ '^l^l’  • 

Using  these  expressions  In  (8.11),  one  can  carry  out  the  differentiations  with  respect 
to  AC  “ Aj^  and  y « p^,  and  one  obtains  (8.15)  for  1 = 1.  The  same  result  is 
obtained  for  1 » n. 

Finally  if  j (2  j ^ n - 2)  is  such  that  sgn9*'(3*  j^)  = -3qn6*'(s*)  = 1,  say, 
(hence  fc'*  has  an  inflection  point),  then  by  (8.8111),  ®g(s)  also  changes  sign  in 
(s^  j,Sj),  and  (8.12)  is  replaced  by 


(8. 17) 


B.(a._i,a^,A,,p^) 


"i-1  ®j 
"i-i  ^ 


(r , cos)  - b , 
1 3 


(k  ^sin)  - d, 

j 3 


where  r.(B,)  = 0,  “j  ^ ^j'  ^ 'A'°  differentiate  the  improper  integrals  one 

replaces  the  A^y  functions  analogous  to  (8.16),  then  (8.11)  for  Ac  * A^ 

and  y = p^  again  yields  (8.15)  for  i = j.  It  should  be  observed  that  depends 

on  a,  ,,ci.,A,,p.,  but  3F  /Sy  and  3G  /3y  do  not  contain  terms  3B./3y.  We  have 


j 


j 


now  established  (8.15)  for  1 « l,...,n. 

We  next  choose  (i  » 1, . . . ,n  - 1)  for  y in  (8.11)  and  obtain 

(Pj^coBo^  + Oj^slnOj^)Kj^^(a^)  - (p j^^^coso^^  + o^^^sina^lx^^^  (a^) 
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I 


or,  usinq  (8.13)  and  (8,15):  k . (ci , ) » k (a  ) , i.e. 

1 1 1+1  i 

(8-18)  0'(s,  - 0)  = 0'(s.  +0),  i = - 1 . 

(11  (XI 

Furthormore,  by  (8.4),  O*(0)  « 0*(s  ) * 0.  Thus  we  have  shown  that  if  (i  is  a 

a an 

critical  i>oint  of  U*(a)  then  0 satisfies  (8.1)  and  (8.2),  hence  0 « Q*,  a = a*. 

a a 

That  conversely  U*(u*)  = U^(/'>)  is  a critical  value  of  U*  follows  immediately  from 
the  fact  that  ^ stationary  value  of  U^.  Proposition  8.2  is  proved. 

The  stability  function  U*  attains  a minimum  in  the  compact  set  N(a*),  say 

U* (a  . ) = min  U* (a)  . 

»./  *v 

(If  N(a*) 

If  a . is  a critical  rxiint  of  U*  (i.e.  a . is  in  the  interior  of  N(a*))  then, 
min  min 

by  the  preceding  proposition,  a . * a*  and  K*  minimizes  the  rxitential  energy  IK 

min  0 

among  all  E with  a f N(a*).  The  theorem  below  will  show  that  in  this  case  E* 
u 

minimizes  n among  all  the  {p^»p,r...#p  )-  interpolants  sufficiently  close  to  A’*, 

0 0 1 n 

hence  that  E*  is  stable.  On  the  other  hand,  if  U*(a*)  is  not  a local  minimum  of 

U*  then  there  are  interpolants  arbitrarily  close  to  K*  for  which 

(£.’  ) = U*(a)  < U*(a*)  = (£’*),  hence  E*  is  unstable.  Thus,  we  arrive  at  the 

0 a 0 

following  effective  stability  criterion: 

Theorem.  SupF>ose  the  indecomposable  extremal  intorpolant  E*  ^ B'  . has  only  proper 

a 

subarcs  /-.’J.  Then  E*  is  stable  if  and  only  if  the  stability  function  U*  has  a 
local  minimum  at  a*. 

Proof.  The  proof  depends  critically  on  the  following  result  which  we  formulate  as 
a lemma. 

Lemma.  There  exists  a neighborhood  N_(0*)  C N(0*)  such  that  (C)  > ) for 

o o — v)  a 

each  C with  normal  representation  0 ^ N^(0*)-  Here  a ® ,}f 

o 1 n-i 

“i  ” 9 ) • 

Proof  of  Lemma.  Since  each  internal  (terminal)  arc  of  E between  consecutive 
a 

interpolation  nodes,  if  considered  as  a 2-point  extremal  interpolant  with  two  (one) 

angle  constraints,  is  stable  it  is  true  that  U_ (C)  > U_(£  ) for  C sufficiently 

0 “ 0 Ql 
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clo.so  to  ,a  fixofi.  'Hio  lemma  assert?;  that  this  inequality  holdji  in  a neiqhborhood 

(t 

that  is  independent  of 

We  may  assume  0 in  the  form  (2.5); 

(H.1*^J)  0 = 0 + tn  + 

ri 

with  h c V (0  ) , ri  (s.  (0  ) ) =0  {i  = 1,  . . . ,n  - 1)  , d°(0,n)  ^ 1,  d°(0,C)  ' 1 . We  also 
0 u 1 a — — 


s.  {0  ) 

I ft 

truiv  assume  d.  * f sinO  / 0 (the  integral  is  independent  of  a),  otherwise  d, 

^ s.  , (0  ) '*  ^ 

1-1  a 

should  be  replaced  by  b^.  Then  by  (2.10),  (1.5) 


/ 0'^  = f + c Q(f)  ,ri)  + R(e), 


(H.20) 


Q(f)^^,n)  = / (O'  - 2 ^ ^ J-  ^ 


s . (0  ) 

1 a 

^ (cosO  )n  / 

n * 


S.  (0  ) 
1 a 


a’n  . 


i=l  s , , (0  ) 
1-1  a 


s.  , (0  ) 
1-1  a 


wiiorf;  K(e)/8.  '0  as  t 0,  uniformly  for  0 e N(0*),  a f N(a*).  The  mappings 

rt  >>•  s (0  ) (i  = 1 n)  , a t'  0 , from  N(a*)  to  IB  , N(0*),  respectively,  are 

1 a ft 

continuous,  and  so  is  the  mapping 


Since  > 0 for  each  a f N(a*)  it  follows  that  = ^*^'^acN(a*)^a  ^ ^ and,  by 

(8.19),  / 0*^  ^ / 0^^^  + ~ t^9(0^,ri)  for  all  sufficiently  small  e,  say  \c\  ^ e^. 

We  can  now  choose  the  neighborhood  N^(0*)  ^ N(0*)  so  that  0 e N^(0*)  may  be 

represented  in  the  form  (8.19)  with  lc|  ^ E^.  Then  / 0’^  ^ / 0^^^,  which  proves 


the  lemma. 

Proof  of  the  Theorem.  We  need  to  prove  only  the  sufficiency  of  the  condition.  Thus, 
we  assume  there  exists  E^  > 0 such  that  U*(fi*)  ^ U*((i)  for  ja  - a*|  ^ e^^.  If 
the  neiqhborhood  C N^^iO*)  is  sufficiently  small  then  |a  - fi*l^Ej  for  each 

0 c Nj(0*),  ft  = (0(3^(0)).  Using  the  Lemma,  we  have 


2 *^2^2 
/ 0*’  . U»(a*)  < U*(a)  - / O'  ^ / 0'  , 

0 0 “ 0 
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^ local  miiiimiun. 

Wc  prf'ifnt  two  cx.JtnplcH,  which  illustrate  tho  effectiveness  of  the  pro{>o:.itions 
in  this  section. 

1_.  Suppose  w<?  have  t he  con f i qurat ion  ^1^0 ' ' 1^2  ^ where*  p^^  * (0,0), 

P^  ■ n,0),  p^  - (J,rj),  Without  loss  we  may  assume  d ^ 1.  It  is  easy  to  see  that, 

lot  each  d,  there  is  an  extremal  interix>lant  A’*,  which  makes  the  anqle  u*  with 

the  v<.*ctoi  at  , where  a*  varies  from  tt/4  to  0 as  d varies  from 

1 t/>  Mer<!  the  stability  function  U*  is  a function  of  a sinqle  variable  u, 

wliich  has  been  computed  by  br.  1).  Pence.  It  is  found  that  U*(a*)  is  a local  minimum 

Jor  each  d.  hy  the  'nieorem,  the  above  A’*  is  a stai)le  extremal. 

lixample  2.  SupiKj  e tlie  conf i<juration  to  be  interix^lated  is  {p,#P^»P^»p  ) with 
' — - 12  3 4 

pj  = (a,t»)»  (1/b),  p^  (0,1),  * (0,a),  where  -“  < a < 1 , This 

conf iqurat  ion  with  a - O.h  was  mentioned  first  in  the  note  |5)  as  an  example  for 

which  there  is  no  intt^rpol.it  iinj  elastica,  and  this  claim  was,  without  examination, 

ti  fieatc'd  It,  many  siil>se<pjent  pul>l  ications.  However,  there  are  interfiolatinq  Qlastica, 

fur  c*arh  a,  in  j^articular  tti«*re  is  one  which  is  symmetric  with  respect  to  the 

symmetry  axis  of  t fie  c’onf  i qurat ion . Tliis  can  be  seen  as  follows.  Let  C„  be  the 

n 

,ynufM!tric  interjx>lant  tjf  ^^^1 '^*2'*’ j'^^4  * wfiich  is  uniquely  defined  by  the  followinq 

cf>ndit.ions.  lias  continuous  slope;  the  arcs  3A  interpolation 

nodes  are  simple  elastica;  A'.,  and  A’_„  have  curvature  0 at  p,  and  p. 

Ip  3p  14 

respectively;  tfie  tanqent  vector  alonq  A'^^  turns  throuqh  the  anqle  fi.  Clearly, 
for  p.  ■ 71,  tlie  curvature  at  p^  jijm[)H  from  0 to  a negative  value;  and  for  some 
p.  ■ 11  the  curvature  at  jumps  from  a negative  value  to  0.  Therefore  there  is 

some  value  (it  if»  unique)  P,^#  0 < < ti,  such  that  ,*C  has  continuous  curvature 

at  p,^  (tfius  also  at  p^,  in  fact  everywhere),  and  tliis  is  an  extremal  interirolant 
of  ^*’l ' *^2'**3' *'4  ^ this  way,  for  each  a,  -“*  < a < 1,  a unique  extremal  inter- 
fx>lant.  A*,  is  defined.  W(!  will  see  that  each  of  these  extremals  In  unstable.  T^e 
results  are  ba.sed  on  computations  carried  out  bv  Dr.  D.  Pence. 


If  a ^ a*,  where  a*  ^ -.27,  then  the  terminal  arcs  of  E*  are  improper, 

hence  E*  is  unstable  by  Proposition  8.1.  If  a < a*  then  the  hypotheses  of  the 

above  Theorem  are  satisfied.  Instead  of  the  stability  function  U*  (a)  ==  *)  » 

^ 0 a* 

a = we  use  the  function  of  one  variable  which  is  the  restriction  of  U*(a) 

to  = 3-rt/4  + (i.e.,  we  consider  only  synmetric  perturbations  of  E*)  , 

TTie  computed  results  show  that  a*  is  not  a minimum  point  of  this  function,  hence  E* 
is  not  stable. 

The  question  whether  there  are  stable  extremal  interpolants  for  the  configuration 
(which  would  necessarily  be  nonsymmetric)  remains  open. 


•>.  Stability  of  closed  exirt^nuils  iji_t^^_rj)alat  In^  re^iil^r  polygons. 

An  extremal  interpolant  //  with  normal  representation  s 0(s)  (0  ^ s ^ s) 
is  said  to  1)0  closed  if 


(^>.1)  0(0)  - 0(s),  0'(n)  « 0*(S)  . 

I,et  pQ*  • * * * » I'n-l  •'*  ri’gnlar  n-gon  (n  3) . In  (2,  Sec.  8]  it 

was  shown  that  there  exist  closed  extremals  that  interpolate  the  configuration 


>P  i>P  P„l-  I'l  particular,  there  la  one,  K , which  has  no  inflection 
n~  1 n 0 n 


[joints.  Ixit  s **  ^ ^ n,  be  its  normal  representation.  Its  total 

O O 

variation  Va(0  ) is  minimal,  Va(0  ) - 2ii.  We  prove 
n n 

F'roiJObi t itjti  0.1.  The  closed  extrc*mal  A'^  (n  ^ 3)  is  stable, 

o 

Proof.  The  course  hi  consists  of  n conqruent  arcs,  each  of  lenqth  1,  and  the 
— n 

o 

increment  of  angle  alonq  each  arc  is  2Ti/n.  We  write  0 for  its  normal  representation 

o o 

ami  define  0 (s  t n)  = 0(s).  Then 

(0.2)  0(s)  » 0(s  - 1)  + 2Ti/n  . 


We  assume;  p^  = (0,0),  and  set  ^ 

d » d > 0.  Then 


(b  ,d  > (k  1, . . . ,n)  with  b « 0, 
k k 1 


f*'  - f''  ° 

(9.  3)  b * / cosd  * -d  sin(k  - l)27F/n,  d,  = / sinG  - d cos(k  - l)2n/n  . 

k-i  k-i 


BecauKu  of  symmetry  wc  have 
(9.4) 

Also, 


0(0)  “ Ji/2  ' 7i/n,  0(1/2)  = 11/2  . 


(9.0) 


(2A)‘^^  - 2 / 


^•^(s)  ■=  AsinO(s),  0 ^ a ^ n 


-1/2  , ,,t,l/2^  , cos'/^'u  du 

os  u clu,  (2A)  d “ 2 J 

0 


Thu  quadratic  form  (3.5)  becomes  in  this  case 

g(0,n) 


I I ln'^(t  + k)dt  - ^'^(t  + k)n^(t  + k)ldt 
k=0  0 


n-1 


- 2 I (A/d  )(/  n(t  * k)co.se(t  + k)dt)‘ 


k=0 
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By  (9.2)  and  (9.5) 


° ° 

( 1/d  ) / n(t  + )c)cos9(t  + )i)dt  = (1/d)  / n(t  + l<) COS0  (t) dt, 
“ 0 0 


(9.6)  Q(e 


n-1  ( 1 

.n)=  I / 

)i=0  ^0 


[n  ' ^ (t  + lc)dt  - ^ {t)n^  (t  + )c) ) dt  - (2X/d)  / n (t  + )c)cos0  (t)dt 


This  form  is  to  be  minimized  on  the  space  (3.4); 


V^(0)  = { 


o o o 

9 ' ''^1  / (b  COS0  + <lsin0)  = 0,  Ic  = l,...,n}  . 

1,2  )c  k 


Here  ^ denotes  the  2 -space  of  functions  of  period  n.  Using  (9.3),  we  find 


(9.7)  V 


_(0)  = (n  e W ! J n(t  + lt)sin0(t)dt  = 0,  )i  = 0,1, ...,n  - l) 

u 1,2 


Put  n(t  + k)  = nj^(t)()i  = 0,1,. ..,n  - 1).  Clearly  Q(0,n|^)  = 0(0,9^)  and 

o o 

if  Hq  t Vq(0).  If  Q(0,n)  attains  its  infimum  for  9 = 9^,  then  also  for 

n = n = (1/n) (n.  + n,  + ...  + n ,),  and  n has  period  1.  For  n of  period  1 (9.5) 
0 1 n“i 


1 pooo  o ^ ®2 


(l/n)Q(0,n)  = / (n*  - An  sine)  - (2nx/d) (J  ncos0) 

0 0 


and  (9.7)  requires  J gsinO  = 0.  Thus,  g must  change  sign  in  (0,1)  and  we 
0 


^22^  2/  ^ 22 

(9.9)  f n’  /n  ^ f sin2irt/dt)  jj  {sin2iTt)  = 4it 

0 0 / 0 

From  (9.5)  we  have  the  estimates 

“ 1/2  -1/2  ' 

< (2iT/n)cos 

(9.10) 

d > cosn/n  . 

With  (9.9),  (9.10)  substituted  in  (9.8),  we  find 


(9.11)  (l/n)Q(0,n)  ^411  (1  - (l/n)tan  it/n  - l/2n  cosn/nl  j n^, 

0 

e e 

thus  n Q(0»n)  is  positive  definite  for  n > 3.  By  Proposition  1,  0 is  stable. 
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